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1 Introduction 



Three-nucleon forces (3NFs) are well established in nuclear physics. Although small compared to the 
dominant two-nucleon force (2NF), they are nevertheless needed to gain a quantitative understanding 
of nuclei and nuclear physics. A recent example in this context is the discussion of the 3NF effects in 
proton-deuteron scattering, see e.g. 012]- Other examples are the binding energy difference between 
3 H and 3 He or the saturation properties of nuclear matter. Only in the last decade a theoretical 
tool has become available to systematically analyze few-nucleon forces and consider such fine but 
important aspects as isospin-violation in such forces and in systems made of a few nucleons. This 
tool is the extension and application of chiral perturbation theory to systems with more than one 
nucleon which require an additional non-perturbative resummation to deal with the shallow nuclear 
bound states and large S-wave scattering lengths. While 3NFs in the isospin limit have been analyzed 
in some detail, see e.g. (HJ |1J EJ E], the question of isospin-violation in the 3NF has not yet been 
addressed in this framework. The work reported here is intended to fill this gap. 

Without further ado, let us address the issues considered here. First, we generalize the classification 
of the isospin dependence of two-nucleon forces due to Henley and Miller [7j to the case of A nucleons 
(A > 3), with particular emphasis on the three-nucleon system, see Section^ This is essentially a 
quantum-mechanical exercise and does not reveal any of the underlying dynamics. The keywords here 
are isospin mixing, charge-independence (breaking) and charge-symmetry (breaking). We stress that 
while such language, which precedes QCD and originates from Heisenberg's definition of isospin to 
account for the almost degeneracy of the proton and the neutron combined with almost equality of 
their strong forces, is useful to categorize few-nucleon forces, in QCD the underlying broken symmetry 
is isospin of the light up and down quarks. This symmetry is broken in pure QCD by the light quark 
mass difference and further by electromagnetism when external electroweak interactions are considered. 
Thus, in the second part of this work, Section |3J we derive the leading and next-to-leading order 
isospin-violating contributions of the 3NF based on chiral effective field theory (EFT)* 5 . We briefly 
recall the counting rules for the inclusion of strong and electromagnetic isospin violation presented in 
[SI and discuss the pertinent terms of the effective chiral Lagrangian in Section [3.11 We present the 
leading and subleading isospin-breaking contributions to the 3NF in momentum space in Section f3,2| 
followed by a brief estimate of the relative strength of these forces in Section \3. 31 We end with a short 
summary. The appendix contains the coordinate space representation of the isospin- violating 3NF. 

2 General considerations 

This section deals with a novel classification scheme for the isospin dependence of the ^4-nucleon 
forces. To derive this scheme, one does not make any assumption about the dynamics underlying such 
forces but only utilizes their transformation properties under isospin symmetry and charge symmetry 
operations on the level of nucleons. 



5 We eschew here pionless nuclear EFT as it is not the appropriate tool to analyze this particular problem. 
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2.1 Definitions and notation 



The non-relativistic A-nucleon system is described by the Hamilton operator H 

H = H + V W + V m + ... + V m , (2.1) 

where Hq is the nucleoli kinetic energy and V nN represents the potential corresponding to the n- 
nucleon force. The total isospin operator T is given by the sum of the isospin operators t of the 
individual nucleons: 

A 

T = Y,t(a). (2.2) 

a=l 

The total isospin operator T as well as the operators t(i) satisfy the Lie algebra of the SU(2) isospin 
group: 

[T%, Tj] = ieijk Tk , (2.3) 
[U(a), tj(b)] = iS ab € ijk t k (a) , 

with k = 1,2,3. The single-nucleon isospin operators ij(a) can be conveniently represented in 
terms of Pauli matrices t%: 

U(a) = Kiia) . (2.4) 
The charge operator Q is defined for the A-nucleon system as: 

Q = e (j +Ts ) ■ (2 ' 5) 

Since the baryon number and the charge are conserved in nuclear reactions, the operator T3 commutes 
with H even if isospin symmetry is broken. 

Charge symmetry represents invariance under reflection about the 1-2 plane in charge space. The 
charge symmetry operator P cs transforms proton and neutron states into each other and is given by 

m 

A A 

Pes = e i7rTa = J] = J] (ir 2 (a)) . (2.6) 

a=l a=l 

Thus charge symmetry conservation means the equivalence of nn and pp, nnn and ppp, . . ., forces. 
Obviously, charge symmetry is valid if isospin is conserved, i.e. if 

[H, T 2 } = [H, Ti}=0. (2.7) 



2.2 Two nucleons 



The classification of the two-nucleon forces according to their isospin dependence has been worked 
out by Henley and Miller 0. For the sake of completeness, we will briefly remind the reader of this 
classification scheme in what follows. 

The two-nucleon forces fall into four classes: 
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• Class (I) forces V 2N are isospin invariant and can be expressed as: 

y/ N = ai + a 2 t(l)-t(2) , (2.8) 
where ccj are space and spin operators. 

• Class (II) forces, V™ , maintain charge symmetry but break charge independence (i.e. are not 
isospin invariant* 6 ): 

[Vjf, T] + 0, 
[V??, P cs ] = . (2.9) 

The class (II) forces are proportional to the isotensor: 

y/ / N = ar 3 (l)r 3 (2). (2.10) 

It is easy to verify that these forces do not mix isospin in the two-nucleon system and thus 
satisfy, in addition, the following relation: 

[Vjf, T 2 ] = . (2.11) 

• Class (III) forces break charge symmetry but do not lead to isospin mixing in the two-nucleon 
system: 

[Vjfj,T] + 0, 

Wffi, Pes] + 0, 

[V™,T 2 ] = 0. (2.12) 

Such forces have the general structure: 

y / 2 / N / = a(T3(l) + r 3 (2)). (2.13) 

and are symmetric under the interchange of the nucleons 1 and 2. 

• Finally, class (IV) forces break charge symmetry and cause isospin mixing, i. e.: 

Wrv, T] + 0, 
W$, Pes] + , 

[V^, T 2 ] + 0. (2.14) 

They can be expressed as: 

Vw = aifo(l) " r 3 (2)) + a 2 [r(l) x r(2)] 3 . (2.15) 
The operator a 2 has to be odd under a time reversal transformation. 



* Clearly, V u as well as all other considered isospin-violating interactions still commute with the third components 
of the total isospin for the reason explained before. 
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2.3 Three and more nucleons 



Let us now generalize the above treatment to systems with more than two nucleons. Considering the 
commutation relations of the Hamilton operator H with the operators T 2 and P cs , one can distinguish 
between four different cases for isospin-violating forces: the Hamilton operator may commute with 
both T 2 and P cs , with one of those operators or with none.* The problem with such a classification 
scheme is that conservation of T 2 depends on the number of particles. In general, an A-nucleon 
force that commutes with the squared total isospin operator in the ^4-nucleon system, T\ = (t(l) + 
i(2) + . . . + t(A)) 2 , will not commute with the operator T 2 A . For example, all isospin-breaking two- 
nucleon forces, which do not cause isospin mixing in the two-nucleon system, lead to isospin mixing 
in the three-nucleon system. On the other hand, the property of charge symmetry is independent 
on the number of nucleons and suitable for generalization. Thus in systems with more then two 
nucleons it is convenient to distinguish between the following three classes of forces: class (I) isospin 
symmetric forces, class (II) forces, which break isospin but maintain charge symmetry and class (III) 
forces, which break both isospin and charge symmetry. For two nucleons, our class (III) interactions 
obviously include the class (III) and (IV) forces in the classification by Henley and Miller. 

Let us now concentrate on the 3N force and list all possible isospin structures. 



Class (I) forces are isospin scalars and have the structure: 

Vj W = E («f + • <D + 7?* [r(0 x r(j)] • t(*)) , (2.16) 

i+j+k 

where otf ' , 0^ and jj are space and spin operators with the superscripts being the nucleon 
labels. 

Class (II) forces satisfy: 

[V n , T] , (2.17) 
[V n , P cs ] = , 



and can be expressed as 

V n = E (4i*3«*3(j)+/3f[r«XT(i)] 3 r 3 (A;)) . (2.18) 

i^j^k 

The forces in eq. (j2.18|) give rise to isospin mixing in the 3N system except in the following two 
cases: 

„ 123 , „ 213 „ 132 , „ 312 231 , „ 321 

al23 a213 _ o312 ol32 _ o231 o321 ( -, Q \ 

Pn — Pn — Pn — Pn — Pu - PlI ■ l / - Ly J 

Class (III) forces satisfy: 

[Vn, T] + , (2.20) 
[V n , P cs ] + . 



#7 In case of two nucleons, only three of these four cases appear, since there are no forces which commute with P cs and 
do not with the operator T 2 . 
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There are four types of such isospin-breaking forces: 

+^T3(i)r 3 (j)T3(A:)) . (2.21) 



The first three terms in eq. ()2.21J) cause isospin mixing in the 3N system except in the following 
special cases: 

123 , 132 213 , 231 312 , 321 

a in + a IH = a IH + a HI = a IH + a IU , 

ol23 o213 _ o312 ol32 _ #231 #321 

Pin ~ Pin — Pm ~ Pin — Pin ~ Pin > 

l}fi + l}fl = iWi + im = Til! + Till • (2-22) 

The last term in eq. (j2.21|) does not lead to isospin mixing in the 3N system. Notice further that 
the quantities fiijk are time-reversal-odd. 

In what follows, we will perform explicit calculation of the dominant isospin-violating three-nucleon 
forces based on chiral effective field theory. 



3 Isospin— breaking three— nucleon force in chiral effective field the- 
ory 

3.1 Power counting and effective Lagrangian 

Isospin-breaking two-nucleon forces have been extensively studied within effective field theory ap- 
proaches, see e.g. [HJ 1121 El Q3I El; as well as using more phenomenological methods, see 
e.g. |161 117j for some recent references. In the Standard Model, isospin-violating effects have their 
origin in both strong (i.e. due to the different masses of the up and down quarks) and electromagnetic 
interactions (due to different charges of the up and down quarks). The electromagnetic effects can 
be separated into the ones due to soft and hard photons. While effects of hard photons are incorpo- 
rated in effective field theory by inclusion of electromagnetic short distance operators in the effective 
Lagrangian, soft photons have to be taken into account explicitly. 

Consider first isospin breaking in the strong interaction. The QCD quark mass term can be expressed 
as 

££ass =-qMq = ~q(m u + m d )(l-er 3 )q , (3.1) 

where 

raj — m,, 1 

e = — ~ - . (3.2) 

m d + m u 3 

The above numerical estimation is based on the light quark mass values utilizing a modified MS 
subtraction scheme at a renormalization scale of 1 GeV JH]- The isoscalar term in eq. 1)3.1(1 breaks 
chiral but preserves isospin symmetry. It leads to the nonvanishing pion mass, M% = (m u + m,d)B ^ 
0, where B is a low-energy constant (LEC) that describes the strength of the bilinear light quark 
condensates. All chiral-symmetry-breaking interactions in the effective Lagrangian are proportional 
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to positive powers of M 2 . The isovector term (oc T3) in eq. (|3.1j) breaks isospin symmetry and generates 
a series of isospin-breaking effective interactions oc (eM 2 ) ra with n > 1. It therefore appears to be 
natural to count strong isospin violation in terms of eM 2 . However, we note already here that isospin- 
breaking effects are in general much smaller than indicated by the numerical value of e, because the 
relevant scale for the isospin-conserving contributions is the chiral-symmetry-breaking scale A x rather 
than m u + m^. 



Electromagnetic terms in the effective Lagrangian can be generated using the method of external 
sources, see e.g. [191 I2U1 121) for more details. All such terms are proportional to the nucleon charge 
matrix Q = e(l + Ts)/2, where e denotes the electric charge.* 8 More precisely, the vertices which 

contain (do not contain) the photon fields are proportional to Q n (Q 2n ), where n = 1, 2, Since we 

are interested here in nucleon-nucleon scattering in the absence of external fields, so that no photon 
can leave a Feynman diagram, it is convenient to introduce the small parameter e 2 ~ 1/10 for isospin- 
violating effects caused by the electromagnetic interactions. As will be discussed below, three-nucleon 
forces due to virtual photon exchange do not contribute at the leading and subleading orders. We 
will therefore not consider virtual photons in the present work. Notice however that electromagnetic 
effects might be enhanced at low energy due to the long range of the corresponding interaction, see 
[Sj for more details. A systematic study of such effects should therefore be performed in the future. 
For the first step in this direction see |22j . 



In the present study we adopt the same power counting rules for isospin-breaking contributions as 
introduced in [Sj. Specifically, we count 

e~e~-|, (3.3) 

where q ~ M n refers to a generic low-momentum scale and A to the hard scale which enters the 
values of the corresponding low-energy constants. In addition, we keep track of the additional factors 
l/(47r) 2 arising from the photon loops by counting 

e 2 o 4 

— 3- (3 4) 

(4vr) 2 A 4 1 j 

Notice further that contrary to the standard practice in the single-nucleon sector, the nucleon mass 
is considered as a much larger scale compared to the chiral-symmetry-breaking scale for reasons 
explained in [Hj. In this work we adopt the counting rule q/m ~ (q/A) 2 , which has also been used 
in |23| . Counting the nucleon mass in this way ensures that all iterations of the leading-order NN 
potential contribute to the scattering amplitude at leading order (q/A)° and thus have to be resummed. 
The N-nucleon force receives contributions of the order ~ (q/ ' A) v ', where 

v = -4 + 2n 7 + 2N + 2L + ^ V- A< . (3.5) 

i 

Here, L and Vi refer to the number of loops and vertices of type i and n 7 is the number of virtual 
photons. Further, the vertex dimension Aj is given by 

A; = di + -Hi - 2 , (3.6) 

where rii is the number of nucleon field operators and di is the q-power of the vertex, which accounts 
for the number of derivatives and insertions of pion mass, e and e/(47r) according to eqs. 1)3.3(1 . (|3.4|) , 

* 8 Or equivalently, one can use the quark charge matrix e (1/3 + ts)/2. 
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Let us now specify the relevant terms in the effective Lagrangian. In the purely pionic sector, we have 
to take into account the following structures: 



Anr = ^f>X + X+> + C(Ql - Ql) , 



(3.7) 



where F n refers to the pion decay constant and the brackets ( ) denote traces in the flavor space. We 
remark that various LECs appearing in the effective Lagrangian correspond to bare quantities in the 
chiral SU(2) limit. Throughout this manuscript we will not specify this explicitly and use physical 
values for the LECs to express our results for the 3NF. Mass and coupling constant renormalization 
is detailed e.g. in refs. (SUES]- Further, 



U, 



X = IBM 



X± = U ^XU ] ± ux*u , 



1 



Q± = ^(u^Qu ± uQvt) 



(3.8) 



The unitary 2x2 matrix U in the flavor space collects the pion fields. In the a-model gauge, it takes 
the form 



TV + IT • 7T 



(3.9) 



The pion mass resulting from eq. (|3.7|) is given by 

B(m u + m d ) , 



B{m u + m d ) + j^e C . 



(3.10) 



The experimentally known pion mass difference M w ± — M n o = 4.6 MeV allows to fix the value of the 
LEC C, C = 5.9 • 1(T 5 GeV 4 . Notice that the natural scale for this LEC is F 2 A 2 /(4it) 2 ~ 3 • 1(T 5 
GeV 4 if one adopts A ~ M p . 



Utilizing the heavy baryon framework, the relevant structures in the single-nucleon Lagrangian are 
[2*6*] (for a more detailed discussion see e.g the review [2*7]^: 



J2 n 



N 



iv • D + gA S ■ u 

+ c 1 <x + > + |<u-u> + ^[S", S"][ Ufl , u v ]+c, X+ 
+ hF 2 (Q 2 + - Q 2 _) + f 2 F 2 (Q + )Q + + f 3 F 2 (Q + ) 2 ]N v 



(3.11) 



where N v refers to the field operator of a nucleon moving with the velocity v^, c\ 345, /i,2,3 are the 
strong and the electromagnetic LECs, respectively, and 



D„ 



^[u\ df.u], X+ = X+~\{x+) , S^ = ]-i^ b a^ v v l - 



(3.12) 

Keeping the terms with at most two pion fields and switching to the nucleon rest-frame system, the 
Lagrangian density in eq. (|3.11j) can be expressed in a more convenient form:* 9 

1 



N 



Art 



ido — Am - 



9A_ 
2F n 



TO ■ V-7T 



4*2 



T • (n X 7f) 



^M 2 tv 2 + J|(*Vr • d"7v) - ^e ijk e abc a t T a (Vj ir b )(V k vr c ) - ^eM 2 {ir ■ t)tt 3 



+h e 2 (ir 2 3 - iv 2 ) + -h e z ((7v ■ t)tt 3 - tt z t 3 ) 



N 



(3.13) 



# 9 Notice that only terms with three and more pion fields depend on the specific parametrization of the matrix U. 
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Notice that at the order we are working, there is no need to distinguish between M^o and M n ± in 
eq. p.l3|) . We have therefore used the same symbol M n for both charged and neutral pion masses. 
The nucleon mass shift Am in the above equation is given by 

Am = -Ac x Ml - \Fle 2 (2h + h + 2/ 3 ) - ^n(4c 5 eM* + f 2 e 2 F 2 ) . (3.14) 

The isospin invariant shift given by the first two terms in eq. (|3,14|) is of no importance and can 
be absorbed by a redefinition of the bare nucleon mass. The proton-to-neutron mass difference 
bra = m p — m n fixes the values of the LECs C5 and f 2 through 

(5m) stT - = (m p - m n ) str - = -4c 5 eM 2 = (-2.05 ± 0.3) MeV , 

(5m) cm - = (m p - m n ) cm - = -f 2 e 2 F 2 = (0.7 ± 0.3) MeV . (3.15) 

These values are taken from [25j. The electromagnetic shift is based on an evaluation of the Cottingham 
sum rule. In principle, this contribution could also be evaluated in chiral perturbation theory including 
virtual photons. While the formalism exists (see e.g |19 [ I20[ l2T]). there are still some subtleties to be 
addressed |29| . Therefore, we consider the electromagnetic mass shifts for the ground state baryon 
octet collected in [2Hj the best values available. Notice that according to the counting rules ()3.3)) and 
(|3.4[) . the strong and electromagnetic shifts in eq. (|3.15j) are effects of order and q A , respectively. 
While the constants C5 and f 2 can be fixed from eq. 1)3.15)1 . the value of the LEC /1, which contribute 
to isospin-violating ttttNN vertex, see eq. ()3.13)) . is unknown. This term plays an important role 
in the analysis of isospin violation in pion-nucleon scattering and the evaluation of the ground state 
characteristics of pionic hydrogen, see |2S] and )32j . respectively. In the two-nucleon sector, it only 
leads to an isospin-invariant contribution to the TPEP at NNLO, which has so far not been considered 
(it can be absorbed in the normalization of the term ~ c\). On the contrary, the resulting contribution 
to the 3NF is isospin-breaking. It, however, does not violate charge symmetry and, therefore, does e.g. 
not contribute to the binding-energy difference of 3 H and 3 He. We further stress that the /1— term 
has not been included in the Lagrangian used in jHUl El EH] since another power counting for the 
electromagnetic effects was employed (see the discussion in 9 ). 

In the few-nucleon sector we only need the following isospin invariant structures 

Cnn = ~C s (N v N v ){N v N v )+2C T (N v S fl N v )(N v S>*N v ) 

~D(N V N V )(N V S ■ uN v ) - ^E(N v N v )(N v tN v ) ■ (N v tN v ) . (3.16) 

where Cst, D and E are the corresponding low-energy constants. The Lagrangian density ()3.16)1 
gives rise to the following relevant terms in the nucleon rest-frame system: 

C NN = -^C s (N j N)(N^N) - ^C T {N^aN)(N^aN) 

- 7 ^-(iV t Af)(iV t o : TiV) • V7r — —E (N^N)(N^tN) ■ (N^tN) . (3.17) 
AF W 2 

3.2 Three— nucleon force in momentum space 

We are now in the position to discuss the leading and subleading isospin-breaking contributions to 
the 3NF # 10 . For the sake of completeness, we will briefly remind the reader of the structure of 

* 10 After submission of our manuscript, related work on charge symmetry breaking in the 3N system by Friar, Payne 
and van Kolck |33| appeared. 
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the isospin-conserving 3NF. The leading 3NF contribution of the order (q/A) 2 represented by the 
graphs in Fig. ^ is we h known to vanish. More precisely, the first two graphs (a) and (b) in this 
figure vanish in the static limit if one adopts an energy-independent formalism such as the method 
of unitary transformation [31] . Alternatively, one can use old-fashioned perturbation theory to derive 
a corresponding energy-dependent 3NF potential. The latter is known to cancel against the recoil 
corrections to the 2N potential being iterated in the scattering equation [SSI 15 ■ ^ should be understood 
that the first two diagrams shown in Fig.^only specify the topology and do not correspond to Feynman 
graphs. Clearly, the corresponding contributions to the 3NF do not include the pieces generated by the 
iteration of the 2NF. We remind the reader that the operators associated with these diagrams depend 
on the scheme and on the definition of the potential. In the method of unitary transformation, these 
graphs subsume both irreducible and reducible time-ordered topologies. The reducible diagrams 
do, however, not contain anomalously small energy denominators, which correspond to the purely 
two-nucleon intermediate states in old-fashioned perturbation theory. The last diagram in Fig. ^ is 
suppressed by a factor of q/m due to the time derivative entering the Weinberg-Tomozawa vertex in 
eq. (|3~T3l . 




(a) (b) (c) 

Figure 1: Leading contribution to the 3NF at the order (q/A) 2 which vanish, as 
discussed in the text. Solid and dashed lines are nucleons and pions, respectively. 
Heavy dots denote the leading-order vertices with A, = 0. 



The first nonvanishing 3NFs arise at order (q/A) 3 from the diagrams shown Fig.[2]with one subleading 
vertex of dimension Aj = 1. The contribution from the first graph in Fig. [2] is also incorporated in 
various phenomenological models like e.g. 
related discussion): 

2 



V 3N 



9A 



the TM99 3NF [H] and given by g] (see also for a 
{Si ■ qi)(aj ■ qj) 



2FJ (e + M 2 )(q, 2 + M 2 ) 



ijk ' 



J3 



(3.18) 



where qt = pi' — pj; pi (pi ') are initial (final) momenta of the nucleon i and 



f: 



a/3 



5 af3 



Ac!M 2 , 2c 3 



F 2 



Here and in what follows, we use the usual notation to express the nuclear force: the quantity is 
an operator with respect to spin and isospin quantum numbers and a matrix element with respect to 
momentum quantum numbers. Notice also that we have changed the notation of section E] and write 
the nucleon labels as subscripts of the spin and isospin matrices (i.e. use Tj and <Tj instead of r(i) and 
a(i)), while the superscripts denote corresponding vector indices. Further, 



ci 



Cl + 



e 2 F 2 h 
2M 2 



(3.19) 
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Note that this renormalization of the sigma-term related LEC c\ by the electromagnetic LEC f\ was 
already discussed in the analysis of pion-nucleon scattering [^S] . Clearly, this electromagnetic shift of 
the LEC c\ represents a higher-order effect and only needs to be taken into account at order (q/A) 5 
and higher. The remaining contributions from graphs (b) and (c) in [51 are given by jH] 



3N 



vr 



9A n Oi-qi 



(3.20) 



(a) 





Figure 2: Subleading contribution to the 3NF at the order (q/A) s . Solid rectangles 
refer to vertices with Aj = 1. For remaining notation see Fig. ^ 



First isospin-conserving corrections to the 3NF arise at order (q/A) 4 , where one has to consider tree 
diagrams with one vertex of the dimension Aj = 2 as well various one-loop diagrams with the leading 
vertices. Derivation of these corrections to the 3NF will be published elsewhere. The main focus 
of the present work is related to isospin-breaking corrections which first appear at the same order 
(q/A) 4 and are given by the graphs in Fig. |31 The first two diagrams (a) and (b) and the last one 
(d) are due to strong nucleon mass shift and of the order e(q/A) 3 ~ (q/A) 4 . It should be understood 
that the proton-to-neutron mass difference has to be taken into account not only for intermediate 
but also for incoming and outgoing nucleon states. The corresponding corrections to the two-nucleon 
force have been recently studied in |14M15| . In what follows, we will not separate the electromagnetic 
and strong shifts in the nucleon mass and express the result in terms of the proton-to-neutron mass 
difference 5m = m p — m n . We use the method of unitary transformation as detailed in 37_ to calculate 
the relevant 3NF contributions. Utilizing the notation of this reference, the corresponding two-pion 
exchange potential can be written as: 



V: 



2rr 



V 



V 

2 



A 1 



Hq — E, 



A 1 



If I 

A 1 



Hq — Efi)(H$ — E v 



-Hi 



Hq — E, 



2 



(Hq — Ef!)(HQ — E Tj 



-H, 



Hq — E v i)(H 
A 1 



E 



-Hi rim 



-H 



n) 
A 2 



H v 



(H - 
A 1 



Er 



-Hi 



(Hq — E v ) (Hq — E v ) (Hq — E ri 



-H x 



7] + h. c. 



(3.21) 



Here rj, rj and fj denote the projectors on the purely nucleonic subspace of the Fock space, while A* 
refers to the projector on the states with i pions. Further, H\ is the leading irNN vertex corresponding 
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to the third term in the first line of eq. Q3.13[) . Hq denotes the free Hamilton operator for pions and 
nucleons corresponding to the density 

n = \k 2 + ^(Vtt) 2 + \mI-k 1 + ~7VW 3 jV , (3.22) 

and E^, and E^ refer to the energy of the nucleons in the states rj, rj and fj, respectively. Notice 
that the first three terms in eq. ()3.21j) subsume the contributions of the reducible graphs while the last 
term refers to the irreducible topology. Neglecting the proton-to-neutron mass difference in eq. ()3.22j) 
one recovers the isospin symmetric result of |37j : 



T,^it~y)Hi f}H\ — H\ + -Hi — Hi fjHi— — kHi — Hi — Hi ■ H i N9 i?i 

2 [UJ) Z UJ I ijj (uj) z LO UJ\ + LOl (iJj z 



ri, (3.23) 



where u denotes the pionic free energy. We remark that eq. (|3.21j) can also be used to calculate 
relativistic 1 / m-corrections to the two-pion exchange potential if one keeps the nucleon kinetic energy 
term in eq. (|3.22|) . An additional unitary transformation should, however, be performed in order to 
end up with the potential used in [S]. 




(a) (b) (c) (d) 

Figure 3: Leading isospin-violating contribution to the 3NF at the order (q/A) . 
Crossed circles refer to isospin-breaking vertices with Aj = 2. For remaining nota- 
tion see Fig. ^ 



Explicit evaluation of the 3NF using eq. ()3.21j) leads to the following result: 



V 3N 
v 2tt 



25m 



9A 



2FJ (% 2 + M2) 2 (^ 2 + M2) 
+ Qi ■ Qj [{ri ■ r fc )r| - (n ■ Tj)rl] X . 



Qjl ■ °k [Ti x Tj\ 



(3.24) 



Notice that we have expanded the energy denominators in powers of 5m in eq. (|3.21|) and kept only 
the linear terms. Similarly to the case of the two-nucleon potential ^1], the resulting 3NF is entirely 
due to irreducible diagrams. As a cross-check of our approach, we have also calculated the two-pion 
exchange 2NF corresponding to eq. ()3.21|) and recovered the results of |14j . 

The contribution of the one-pion exchange diagram (b) in Fig- El is given by the operators 



Vin = rf 



A 1 



— H i 

2 (Hq — E v )(Hq — Ef f 



1 



-H 1 fjH 2 + -H 



Ho 



2 (Hq — E v ) (Hq — E, t 



-Hi 



rj + h. c. (3.25) 
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where H 2 corresponds to the first two terms in eq. (|3.17j) . Similarly to the previously considered case, 
we recover the result of [37] in the limit 5m — * 0: 



Vi 



1 A 1 1 A 1 A 1 A 1 

-Hx—H^K?, - -HifiHx—Hx + U x —R 2 —H x 



We find the following expression for the isospin-breaking one-pion exchange 3NF 



V x f= 2 5mC T 

i+j+k 



9A 



[Tk x Ti\ [<Tj x a k \ ■ qi . 



2FJ (% 2 + M 2 ) 2 

Notice that can be rewritten in an equivalent form making use of the relation 
[r k x Ti} 3 [3j x a k ] ■ q} = ({n ■ Tj)t| - (r, • r k )Tj) {Bj ■ $) , 



(3.26) 



(3.27) 



(3.28) 



which holds true when the corresponding operators act on antisymmetrized states with respect to j 
and k. 

The diagram (c) in Fig. |3]is due to the cs-term in eq. (|3.13|) and of the order e{q/A) 3 ~ (g/A) 4 as 
well. Denoting the interaction oc C5 by H3, the contribution of this graph is given by 



V: 



2tt = V 



X 2 



A 1 A 1 A 1 
Hi — #3 — Hi + Hi — Hi- ■ - 

U UJ UJ [UJl + U)2, 



H 3 + H 3 



X 2 



(3.29) 



' (u>i + 0J2) w 

Alternatively, one can use the Feynman graph technique to evaluate the corresponding 3NF. We find 



Vn 



3N 



2tt 



E 

i^j^k 



{5m 



istr. 



9A 



{Si • qi)(vj ■ qj) 



4F2 \2Fj {qf + M 2 ){qj 2 + M%) 



{Ti • T k )T. 



(3.30) 



Notice that all leading (i.e. ~ {q/A) 4 ) isospin-violating 3NFs given by eqs. (j3.24|) . (|3.27|) and (|3.3()|) 
are charge-symmetry-breaking, i.e. of class (III) in the notation of section 12.31 We further point 
out that although the M w ± 7^ M w o -corrections to the graphs in Fig. ^ given by the first three graphs 
(a), (b) and (c) in the later Fig.|f)]are formally also of the order {q/A) 4 , they lead to 1/m-suppressed 
contributions to the 3NF for the same reason as do the corresponding isospin-conserving terms. 



The contribution of the last diagram (d) in Fig. |21 is given by 

A 1 „ A 2 



2 ' 



{Hq — E v ) {H — Efj) 



-H, 



A 1 



+ #1 



A 1 



{Hq — E Tj 



-HI 



WT 



A 1 



{Hq — Erj 



Hi 



{Hq — E v ) {H — E v ) 
rj + h. c. , 



H, 



(3.31) 



where H^ T refers to the Weinberg-Tomozawa vertex. Explicit evaluation of this graph can be per- 
formed expanding the above expression in powers of 5m and keeping the terms oc 5m. Alternatively, 
one can use the Feynman graph technique. In that case one should use for the energy transfer of the 
nucleon i: q® = {p'j} — = Am + 0(m _1 ), where Am denotes the nucleon mass difference in the 
final and initial state. We find: 



T/ 3N 

V 2tv 



E 



5m ( g A 



{ffi ■ qj){ffj ■ qj) 
4F 2 V 2F„ J («f 4 2 + M%)(qj 2 + M 2 ) 



[{Ti • T k )rf - {Ti ■ Tj)4] . 



(3.32) 
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- -0- 



(a) 




(c) 




(c) 



(f) 



Figure 4: Subleading isospin- violating contribution to the 3NF at the order (q/A) 5 . 
Crossed rectangles refer to isospin-breaking vertices with Aj = 3. For remaining 
notation see Fig. ^ 



The first corrections to the leading isospin-breaking 3NFs arise from the diagrams (a) (b) and (e) in 
Fig. Eland are of the order (e/4-7r) 2 g/A ~ (c//A) 5 . Notice that the contributions of the graphs (c), (d) 
and (f) in this figure are already included in eqs. (|3.24j) . I|3.27|) and (|3.32|) . The first two graphs in 
Fig. 0] represent isospin-violating corrections to the graphs (a) and (b) in Fig. [21 due to the pion mass 
difference and lead to 

T/3N V- KM* ( 9A ^ ft ' IM? ' M f-3_3 [ , 2c 3 -j x 



i^j^k 



+ ^rf [tj x T k f [qt x qj] ■ a k 



< = " E S ^£* D m^ ^Vr®, (3-33) 



i^j^k 

where we have defined 



{e + M%y 

5Ml = Ml ± - M%> . (3.34) 

Notice that at this order (i.e. at (g/A) 5 ) one has to distinguish between the charged and neutral pion 
masses in the pion propagators in eqs. (|3.18|) and l)3.2Uj) . Isospin-violating corrections in eq. ()3.33|) 
are consistent with taking M w ± in the pion propagators in eqs. (|3.18|) and (|3.20|) . The contribution 
of the diagram (e) can be obtained from eq. ()3.29|) : 

+ f ie 2 rfrf). (3.35) 



The 3NFs resulting from M w ± ^ M n o in graphs (a) and (b) of Fig. 0] are charge-symmetry-conserving 
(i.e. class (II)) while the diagram (e) in this figure gives rise to both charge-symmetry-conserving 
(oc f\) and charge-symmetry-breaking (oc (5m) em ') 3NFs. We stress again that the contribution 
~ fl is considered here for the first time. Notice further that the charge-symmetry breaking 3NFs in 
eqs. (|3,30[) . ()3.32j) and ()3.35|) can be combined to: 

(8m) stT - ( g A \ 2 (<?i • qi)(ffj ■ % 



T/3N _ / <M \ K"i • qi)\uj • <4j) r , W 3 , w 31 (a oa\ 



i: _ iL 1/^ \2FJ (& + M2)(tf + M*) 
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The coordinate space representation of the obtained 3NFs is given in appendix [X] Notice that there 
exist further diagrams at this order which, however, lead to vanishing contributions and are not 
considered in the present work. 

Let us now comment on the obtained results. First of all, we notice a (formally) larger relative size of 
the isospin-breaking corrections compared to the two-nucleon sector. Indeed, isospin-breaking 3NFs 
are suppressed by q/A compared to the isospin-conserving 3NFs, while the suppression factor in case 
of the 2NF is (q/A) 2 . Secondly, the leading isospin-breaking corrections to the 2N and 3N forces 
arise from different sources. In particular, the dominant contribution to the 3NF is governed by the 
proton-to-neutron mass difference, which only gives a sub-subleading isospin-breaking correction to 
the 2N force. Further, charge dependence of the pion-nucleon coupling constant does not show up in 
the 3NF at the considered order. Similarly, the leading isospin-breaking 3N contact interaction is of 
the order eM%(q/ A) 3 ~ (q/A) 6 and therefore does not need to be included. Last but not least, we 
notice that the hierarchy of isospin-violating forces observed in the two-nucleon system (i.e. charge- 
independence-breaking forces are stronger than charge-symmetry-breaking forces :9.5) is not valid for 
three-nucleon forces. 

3.3 Estimation of the size of the isospin-breaking 3NFs 

Having derived the dominant isospin-breaking 3NF corrections it would be very interesting to see 
how large the effects actually are. This, however, requires explicit calculations of few-nucleon observ- 
ables, which goes beyond the scope of the present study. Here we restrict ourselves to the following 
very rough estimation. Consider the two-pion-exchange correction given in eq. (|3.24|) . Approxi- 
mating l/(<7i 2 + M 2 ) ~ 1/M 2 we obtain the same spin-space structures as the ones which enter 
the leading isospin conserving 3NF in eq. (|3.18|) . Neglecting the isospin structure one observes that 
the strength of the isospin-breaking terms in eqs. (|3.24|) . (|3.30|) and (|3.32j) reaches few percent of 
the strength of the corresponding isospin-conserving pieces in eq. (|3.18|) . Based on the above esti- 
mates and on the fact that two-pion exchange 3NFs typically contribute several hundreds keV to the 
binding energy of 3 H and 3 He,* 11 one might expect the contribution of the isospin-breaking 3NF in 
eq. (|3.24|) to the 3 He- 3 H binding-energy difference to reach 10 ... 20 keV. On the other hand, the 
relative strength of the formally subleading two-pion exchange terms in eqs. (|3.33[) and H3.35|) reaches 
even 2<5M 2 /M 2 ~ 15%. This surprisingly large size of the subleading isospin-breaking corrections 
compared to the leading ones is due to the LECs C13 4, which enter eq. (|3.33|) and are numerically 
large (the physics behind this enhancement of the LECs is well understood [SHI* 12 )- Notice that a 
similar situation occurs for the isospin-conserving two-pion exchange 2N force, where the numeri- 
cally dominant contributions are provided by subleading terms. One should, however, keep in mind 
that the isospin-breaking 3NFs oc ci 3 4 do not lead to charge-symmetry-breaking and thus do not 
contribute i.e. to the 3 He- 3 H binding-energy difference. The leading charge-symmetry-breaking one- 
pion-exchange 3NF in eq. (|3.27j) is numerically smaller in size than the corresponding subleading 
charge-symmetry-conserving contribution in eq. ()3.33|) as well, although the reason is now completely 

#11 In |38| contributions of various pieces of the Tucson-Melbourne 3NF to the 3 H are considered. While the so-called 
a-term (it corresponds to the ci-term in the chiral 3NF) was found to provide only a tiny contribution, the b- (oc C3) 
and d-terms (oc C4) give about 250 . . . 300 keV each. In the analysis 6| based on chiral EFT, the expectation value of 
the two-pion exchange 3NF for 3 H (with the reduced values of the LECs 03,4) was found to be 390 . . . 730 keV depending 
on the cut-off chosen. 

* 12 In that paper is was shown, that the smallness of the NA mass splitting enhances certain pion-nucleon LECs when 
one integrates out the delta. Furthermore, scalar and vector mesons make large contributions to ci and C4, respectively. 
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different. The 3NF in eq. (|3.27j) is proportional to the LEC Ct which is numerically small j40j.^ 13 
It should be noted in this context that the size of the isospin-breaking 3NFs would be more natural 
if one would treat A-isobar as an explicit degree of freedom. In that case a large portion of the 
subleading 3NFs oc 03,4 and D due to graphs (a) and (b) in Fig. [5] would be promoted to the leading 
order. Note also that such an approach with explicit deltas is much more complicated since one has 
to deal with more structures and also needs e.g. to reanalyze pion-nucleon scattering (for an attempt 
see e.g. |42]). Further, one should keep in mind that the above numerical estimates are very rough. 
In particular, taking into account the neglected isospin structure will change the numbers by several 
times depending on the process considered. Thus, only explicit calculation of various few-nucleon 
observables will provide quantitative insights on the size of the derived 3NFs. 




(a) (b) (c) (d) 

Figure 5: Isospin- violating contribution to the 3NF due to intermediate A- 
excitation (double lines), which are not considered explicitly in the present work. 
The effect of such diagrams is hidden in certain contact operators that originates 
from integrating out the delta in the approach considered here. For remaining no- 
tation see Fig. ^ 

Finally, we point out that there are many 1/m-corrections to the obtained results, some of which 
are depicted in Fig. H3 Since we consider the nucleon mass as a larger scale compared to A, such 
relativistic corrections are irrelevant at the order considered in this work. Notice, however, that if 
one would adopt the counting rule m ~ A, various 1/m-corrections (including the ones due to virtual 
photons) would have to be included at the subleading order (q/A) 5 . Some 3NF diagrams due to virtual 
photon exchange have been considered by Yang and found to provide relatively small contributions of 
the order of ~ 7 keV to the 3 He- 3 H binding-energy difference |431 I44| . Furthermore, we remind the 
reader that the long-range electromagnetic 3NFs might, in principle, give rise to large contributions 
to scattering observables under certain kinematic conditions [S]. 

4 Summary 

Here, we summarize the pertinent results of this investigation. 

i) We have given a classification scheme for j4-nucleon forces according to their isospin dependence. 
In the 3N system, one finds three different classes of forces, according to their transformation 

#13 In EFT without or with perturbative pions, one has Ct = in the limit when both NN S-wave scattering lengths 
go to infinity |41|. 
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(a) 




< -dvrv/wt) 



CWVUVWVD 



(c) 



(d) 



(e) 



Figure 6: Selected 1/m-corrections to the isospin-violating 3NF, which are not 
considered in the present work. Wavy lines refer to photons and unfilled circles 
denote vertices with photons. For remaining notation see Fig. ^ 



properties under isospin and charge-symmetry transformations. 

(ii) We have worked out the leading and subleading isospin-violating 3NFs. The leading contri- 
butions are generated by one- and two-pion exchange diagrams with their strength given by 
the strong neutron-proton mass difference. The subleading corrections are again given by one- 
and two-pion exchange diagrams, driven largely by the charged-to-neutral pion mass differ- 
ence and also by the electromagnetic neutron-proton mass difference and the dimension two 
electromagnetic LEC /i, that plays an important role in the pion-nucleon system. 

(iii) We have estimated the relative strength of the leading and subleading corrections compared to 
the isospin-conserving 3NF at the same order. Isospin-violating 3NFs are expected to provide 
a small but non-negligible contribution to the 3 He- 3 H binding-energy difference. 



In the future, these isospin-breaking forces should be used to analyze three- and four-nucleon systems 
based on chiral EFT, extending e.g. the work presented in [Sj. 
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A Coordinate space representation 



The leading and subleading 3NFs are local and can easily be transformed into coordinate space. We 
first define the following operators: 

d 3 qi d 3 qj - M r {di ■ qi){Bj ■ Bj) 



1 _ f d^d^ ^ iqyr (B l -q l )(B r q J ) 

~ J (2vr)3 (2vr)3 (tf + M^tf + M|) l * ^ J 

= (<?i • V ik )(Bj • V ife )[V ifc x Vjfc] • B k h 2 {r ik ) hi{r jk ) 

= ^2 ~ ~T~ ((^» ' " f jfc) x f id ■ ^fc C 1 + x ^)( 3 + 3x ife + 

+ [(Tj X CTj] • CTfc (1 + - (CTj • f ik ) [f ik X CTj] ■ B k (1 + + Xjk) 

- [di x r jfc ] • B k (°j ■ fjk) (3 + 3x jk + x 2 fe )) 

+ 24tt T^T" ' fifc ) t fifc X ^ ' ° k ^ + Xife ) ~ & X ^ ' ( AA ^ 

2 _ /* d^d^ ^ ^ (B t ■ qj)(Bj ■ Bj) ( _ , 

7 (2vr)3 (2vr)3 (£2 + M 2)2(£2 + M 2) W * ^ 

= ^T~2 r~ ( ' " " ( x + x ik)(3 + 3x jk + x 2 jk ) 

oZTT Vi k T 'jfe V 

+ (<7j • CTj) (1 + Xjk) ~ (Bi ■ f ik ) (Bj ■ f ik ) (1 + X ik )(l + Xjk) 

- (Bi ■ r jk ) (Bj ■ r jk ) (3 + 3x jk + x 2 jk f) 

3 _ f d^d^ i?i . rik i? .. Pjk (Bi-qi)(B 3 -q 3 ) 

^ k J (2vr)3(27T)3 " (£2 + M2)(£2 + M2) 

= -(a, • V ik )(Bj • Vjfc) /ii(r ife ) /ii(r jfc ) 

= T^i T~ (°i " rik)(°j ■ fjk) (1 + x ik ) (1 + Xjk) (A.3) 

167r r 4 r jfc 

Q 4 = f d^d^qj_ ^ {Si-q^ffj-qj) 

U ^ k J (2VT)3 (2VT)3 ^.2 + ^2)2(^.2 + ^) 

= -(<?i • V ik )(Bj ■ Vjk) h 2 (r ik ) hi{rjk) 

~ Xlk ~^2—& " ^fcX^j • (1 + »jfc) ( A - 4 ) 

jr'fc 



32vr 2 r: 



U ^ k J (2tt)3 (2vr)3 (q? + M2)2 ^ ^ 

-(<?i • Vj fc )(a,- • V ik ) h 2 {rik) g(rjk) 

1 e _:Eifc / -> \ 

■ s3 ( r jk) ( (Bi ■ fi k )(Bj ■ fik) (1 + Xjk) - {Bi ■ Bj) J . (A.5) 



8vr rj fe 
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Here is the relative distance between the nucleons i and j, r 
Further, 



ij r ijl> ""ij '13/ ' 13 



r ij/ r ij an d Xij — M^Vij. 



h x {r) 
h 2 {r) 
9(r) 



d 3 q 

(2^)3 (q 2 + M2) ~ 4^' 
d 3 q 



„iq-r 



-M n r 



(2vr) 3 (q 2 + M 2 ) 2 SttM* 



(2vr) 3 



6». 



(A.6) 



The isospin-violating 3NF in eqs. ()3.24j) . 1)3.27)1 . ()3.3U)) ()3.33)) and ()3,35|) can now be expressed in 
terms of the operators Ojj k 5 defined above: 



V 



3N 



2 (H) 6m ~ (5m)Str ' ^ + 2SmC T°U 



3 / 9 A 



9A_ 

2F W 
2 



^ + ^p2(^) Str - °%k ~ 25m C T 5 



2F„ 



25m O 



1 



+ TP [Tj X Tfc] 3 ^ C 4 Oj. fc 



_i_ ^3^3 
+ T i T j 



^5M^0 2 ijk + he 2 Of jk - ^ Cl Ml5MlOf jk - ±- 



D5MlO% k 



(A.7) 



Notice that the expressions for the operators 0\^ k in eqs. I)A.1)I - I)A.5)) are singular at short distance 
and need to be regularized. If one chooses to work with the local regulating functions, the regularized 
expressions can easily be obtained by an appropriate modification of the functions h\(r), hzir) and 
9(r). 
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